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Abstract 

On December 27th 2004, a giant 7-flare from the Soft Gamma-ray Repeater 1806-20 saturated 
many satellite gamma-ray detectors. This event was by more than two orders of magnitude the 
brightest cosmic transient ever observed. If the gamma emission extends up to TeV energies 
with a hard power law energy spectrum, photo-produced muons could be observed in surface 
and underground arrays. Moreover, high-energy neutrinos could have been produced during the 
SGR giant flare if there were substantial baryonic outflow from the magnetar. These high-energy 
neutrinos would have also produced muons in an underground array. AMANDA-II was used 
to search for downgoing muons indicative of high-energy gammas and/or neutrinos. The data 
revealed no significant signal. The upper limit on the gamma flux at 90% CL is dN/dE < 0.05 
(0.5) TeV -1 m~ 2 s _1 for 7 = —1.47 (—2). Similarly, we set limits on the normalization constant 
of the high-energy neutrino emission of 0.4 (6.1) TeV -1 mT 2 s _1 for 7 = —1.47 (—2). 

PACS numbers: 95.55.Vj, 95.55.Ka, 95.85.Pw 
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Introduction. — Soft Gamma-ray Repeaters (SGRs) are X-ray pulsars which have qui- 
escent soft (2-10 keV) periodic X-ray emissions with periods ranging from 5 to 10 s and 
luminosities of the order of 10 33-35 erg/s. They exhibit repetitive bursts lasting ~ 0.1 s 
which reach peak luminosities of ~ 10 41 erg/s in X-rays and 7-rays. There are four known 
SGRs, three in the Milky Way (including SGR 1860-20) and one in the Large Magellanic 
Cloud. Three of the four known SGRs have had hard spectrum (~MeV energy) giant flares 
with luminosities reaching up to ~ 10 47 erg/s. The first of these giant flares (from SGR 0525- 
66 l|) was observed on March 5, 1979 by the Venera 11 and 12 spacecraft. SGR 1900+14 
exhibited a giant flare in 1998 The most recent and brightest flare came from SGR 1806- 
20 on Dec. 27, 2004. This flare lasted about 5 minutes (the duration of the initial spike 
was ~ 0.2 s), had a peak luminosity of ~ 2 • 10 47 erg/s and a total energy emission of 
~ 5 • 10 46 erg 3]. This event was observed by several satellite experiments P, U, 0], al- 
though they saturated during the blast. Recent estimates locate the source at a distance of 
15.1±i;| kpc Q,but this value is still under debate fl. 

The favored "magnetar" model for these objects is a neutron star with a huge magnetic 
field (B ~ 10 15 G). These giant flares can be explained as global crustal fractures due to 
magnetic field rearrangements liberating a high flux of X-rays and 7-rays |8j . 

During the flare, the 7-ray spectrum up to ~ 1 MeV is well described with a blackbody 



spectrum with kT =175 keV 



though the last energy bin may leave room for a power- 



law component. Similar considerations apply to the measurement up to 10 MeV in Ref. 
However, fits to the data favor the presence of a non-thermal component. Fits using a 
blackbody + power law (BB+PL) spectrum show that the PL spectral index starts to 
harden months before the giant flare up to values of about — 1.3 -j — 1.4 |2J below 10 keV. 

The possibility of using underground detectors to observe the muons produced in the 
electromagnetic showers induced by TeV gammas generated in these flares was presented in 
Ref. flfil ]. There have been suggestions of substantial barvonic outflow and the possibility of 
high-energy neutrino production. Radio observations III, |l2| indicate an expanding radio 
source with velocity 0.25-0. 40c. Gelfand et al. argue that a re-brightening of the radio 
emissions ~20 days after the giant flare can be explained if substantial amounts of released 
energy went into a baryonic fireball, and make predictions for TeV neutrino production. Ioka 
et al. ^| also argue that high-energy neutrino production can be related to the fraction of 
burst energy released in the form of baryons. 
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The Dec. 2004 giant flare represents an excellent opportunity to probe the high-energy 
spectrum of these sources by looking for events correlated in time and space with this flare. 
In this paper we present the results of a search for a gamma and/or neutrino signal during 
the SGR 1806-20 giant flare using data from the AMANDA-II detector. The short duration 
of these events and the fact that they come from a point-like source result in ne glig ible 



atmospheric muon and neutrino backgrounds. For the first time in AMANDA (see 15( for 



an example of previous search in other kind of sources) , we used down-going muons to look 
for muon photo-production in the atmosphere, already proposed in Ref. |l6j]. A search for 
coincidences using gravitational waves was presented in Ref. [l^. 

AMANDA (Antarctic Muon And Neutrino Detector Array) is currently running in its 
AMANDA-II configuration of a 3D array of 677 optical modules (OMs)distributed along 
19 strings deployed at depths of 1500-2000 m in the South Pole ice [18]. These 8-inch 
photomultipliers, enclosed in pressure-resistant glass spheres, make it possible to reconstruct 
direction and energy of relativistic muons through timing and intensity of the Cherenkov 
light. The ice layer above the detector reduces the background of atmopheric muons by 
more than 5 orders of magnitude compared to the surface flux. Events are recorded when 
at least 24 OMs register a signal within 2 /xsec. The detector rate on December 27th was 
90 Hz (close to the AMANDA-II average). 

Analysis technique. — For events like the burst of the SGR 1806-20 in Dec. 2004, in which 
most of the energy emitted during the flare is concentrated in a 1 s time scale, the precise 
time and location of the event imply that the background of atmotspheric muons becomes 
negligible thereby allowing a search for TeV 7-rays and downgoing TeV neutrinos. 

The IceCube collaboration follows a policy of blindness in its analysis strategies. By 
studying the expected backgrounds and signals prior to looking at the data, the analysis can 
be designed in an unbiased fashion. In the case of expected small signals, this is particularly 
relevant for having a clear procedure to determine the probability of an event to be produced 
by background. Thus, in this analysis, the determination of the optimum selection criteria is 
done using the simulation of the signal and comparing it with the expected background. The 
procedure, similar to that used in the search for upward moving neutrinos from gamma-ray 
bursts is as follows: 

• the background on-source and off-time is calculated using real data, keeping blind 10 
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min around the burst onset. About one day of off-time data was used to monitor the 
stability of the detector; 

• the signal from the source is simulated in order to estimate the angular resolution and 
the effective of the detector; 

• the appropriate time window is estimated, based on the flare onset times given by 
different X-ray satellites and their counting rates; 

• the optimum search bin size is found by minimizing the Model Discovery Factor 
(MDF) H, 

defined as 

MDF = ^h£M3 (!) 

n s 

where \i is the Poisson mean of the number of signal events which would result in 
rejection of the background hypothesis, at the chosen confidence level CL, in SP% 
of equivalent measurements. SP stands for statistical power and n s is the number 
of signal events predicted by the model. This definition is analogous to the Model 
Rejection Factor (MRF see Ref. Q) 

that is used for setting upper limits. In the case 
of the MDF, the bin size is optimized to maximize the probability of discovery (for 
CL corresponding to five sigma and SP=90%); 

• once the optimum search bin size has been found, the unblinding of the data is done, 
i.e. the events inside the time window and the search bin are counted; 

• this number of events is translated into a flux or a flux limit if no significant excess on 
top of the background is found through the knowledge of the expected signal in the 
detector for the given analysis cuts. 

Although both TeV 7-rays and neutrinos produce muons in the detector array, the optimal 
choice of selection criteria depends on the assumed signal. The analysis was optimized to 
the TeV 7 signal. Any further optimization to the neutrino signal would be more than offset 
by the penalty for an additional trials factor. 

Data and simulation. — In order to have a background estimate for the flare, a time and 
angular window have to be defined around the flare of equatorial coordinates (J2000) right 
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21 



221. The bulk of 



ascension = 18h 08m 39.34s and declination = -20° 24' 39.7" 
the flare energy was concentrated in less than 0.6 s. The Swift-BAT counting rate drops 
by more than 2 orders of magnitudes after 0.6 s from the onset of the burst [9J. Based on 
the observation time for each satellite j^O,^, 23, 24, 2], and accounting for its position, 
the expected signal time in AMANDA was calculated. The spread of the resulting times 
indicate that a safe window is 1.5 s, centered at 21h 30m 26.6s of Dec. 27th, the onset time 
of the flare. 

To evaluate the performance of the detector, simulations were performed for several 
input signals. The CORSIKA-QGSJetOl |26J] and ANIS |23| codes were used to simulate 
the photon (and proton) and neutrino interactions, respectively. The generated energy 
for photons is 10 TeV to 10 5 TeV and for protons is 10 GeV to 10 5 TeV. The muons were 
propagated in the ice with MMC 



28| and the program AMASIM [29] simulated the response 
of the detector. The tracks were reconstructed with the same iterative log-likelihood fitting 
procedure that was applied to the real data. 

The angular resolution for different cuts has been studied using the simulation of down- 
going muons generated by cosmic rays. The angular resolution, defined as the median of 
the angular difference between the true and the reconstructed track, is 3.5°. This value 
was obtained using atmospheric down-going muon high statistics simulations and it was 
checked that this result is robust within 0.1 deg for the expected signals of photons and 
neutrino induced muons. We also considered a variety of spectral indices for the 7-ray 
spectrum assuming values given in Ref. [ljj and we found that the angular resolution is 
almost independent of the spectral index. The effective areas, defined as the equivalent area 
for a perfect detector that is able to detect particles with 100% efficiency, for gammas and 
neutrinos are shown in figure ^ as a function of the energy. 

As noted earlier, the optimum angular window is determined by minimizing the Model 
discovery factor. In figure |2] we show the dependence of the MDF and MRF on the circular 
angular window around the source position. The steps in the MDF curve are due to the 
discreteness of the Poisson distribution. It can also be seen that the MRF minimum interval 
is quite broad, indicating that this variable is not very sensitive to the increase in the number 
of background events with the increase of the angular window. This is due to the small value 
of the background in the allowed time window. 

The search bin size which optimizes the probability of discovery is at 5.8°. With this cut 



8 



E 



10 J 



10 2 



10 



_l I I L. 



4.5 



_l I I L. 



5.5 



6.5 7 7.5 8 

log io E (GeV) 



FIG. 1: Effective area after reconstruction and track quality selection for gammas (solid, black) 
and neutrinos (blue, dashed). 

and a 1.5 s time window, the expected background is 0.06 counts. For events that satisfy 
the detector trigger, we keep almost 80% of the signal in this angular window. 

Systematic uncertainties. — Several sources of systematic uncertainties have been con- 
sidered in this analysis. The uncertainty of the hadronic model calculation is estimated to 
about 15%, mostly due to the unknown fraction of diffractive p production j^J. We have 
compared the muon yields using two different models -CORSIKA/QGSJET and an analytic 
calculation [16]- and found that they agreed to within 5%. The uncertainty in the detector 
efficiency (20%) comes mainly from the overall sensitivity of the OMs and the optical prop- 
erties of the ice. The effect has been estimated simulating different reasonable ice models 
and OM sensitivities. 

The stability of the run was checked in order to exclude possible non-particle events 
induced by detector electronics. These events are identified by a specific method j^l( looking 
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FIG. 2: Model discovery factor (solid, black) and model rejection factor (dashed, red) for an E 1A7 
spectrum. 

for anomalous values in a set of denned variables. A correction is made for the electronics 
dead time (17%, which is a typical value in normal runs). Finally, the simulated and 
measured distributions of an extensive set of variables, like zenith angle, number of hit 
optical modules, and hit times, were compared in the search of possible anomalies. In all 
the cases the agreement was within the systematic errors estimated above. 

Results. — Once the optimum search bin size of 5.8° around the source was determined, 
we unblinded the 1.5 s data around the burst looking for events satisfying the analysis re- 
quirements. No event was observed in the on-source, on-time window. Then, we determined 
the upper limits [32J of the normalization constant A 9Q at a CL of 90% assuming a power-law 
energy spectrum, 

^ < A 90 (E/TeVy (2) 
with a cut-off at 10 5 TeV. These limits are shown in figure H3 together with the sensitivity 
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of the detector. 

To give an idea of the impact of these limits on theoretical estimates such as the 7 flux 
extrapolations presented in Ref. ^(|, for spectral index —1.47 (—2) the limit on the gamma 
flux normalization constant is Aqn = 0.05 (0.5) TeV" 1 m~ 2 s -1 . The calculation of this 

n 

limit using the same energy limits as in Ref. hO would give 3.3 (33) TeV 1 m 2 s 1 , which 
rule out spectral indices 7 ~ —1.5, but not softer (assuming a maximum gamma energy 
of 500 TeV). The effect of the attenuation of the gamma flux by the cosmic microwave 
background and the Galactic interstellar radiation field has been also taken into account 
and has been calculated from the results of Ref. [3^. 

Since the source is above the horizon (hence there is not much column depth for neutrinos 
to interact), the neutrino flux limits are an order of magnitude worse than the TeV 7 limits, 
but can still be used to constrain models. In cases where there is large baryonic outflow, 
high-energy neutrinos are produced and the baryons may make the source partially opaque 
to high-energy photons. Comparing the extrapolations in Ref. [h]], for spectral index —1.47 
(—2) the limit on the flux normalization is A§§ = 0.4 (6.1) TeV -1 m -2 s _1 while the 
model predicts (accounting for oscillations) 1.7 (4.1 x 10~ 4 ) in the same units. We are thus 
able to exclude an extremely hard neutrino spectrum extrapolated from the measured MeV 
photon flux. On the other hand, our limit on the high-energy neutrino fluence is still at 
least one order of magnitude larger than the fluence predicted in Ref. jl^ . 

Conclusions. — In summary, we have searched for TeV gammas and neutrinos associated 
with the Dec. 27th giant flare from SGR 1806-20 . We demonstrate that underground 
neutrino arrays such as AMANDA and IceCube can be used as TeV 7 detectors for transient 
events. An analysis of AMANDA data yields no muons coincident with the flare. We use 
this muon non-observation to place stringent limits on TeV radiation from this source. 
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